Lyapunov exponent for a gas of soft scatterers 
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For a fast particle moving within a two-dimensional array of soft scatterers - centers of weak 
and short-range potential - the dependence of the Lyapunov exponent on the system parameters 
is studied. The use of the linearized equations for variations of the propagation angles and impact 
parameters of consequent collisions reduces the problem to that of calculation of the Lyapunov 
exponent of an ensemble of strongly correlated random matrices with given statistics of matrix el- 
ements. In the simplest approximation this Lyapunov exponent is proportional to the interaction 
strength and inversely proportional to the square root of the interaction range. The model satisfac- 
torily describes the intensity of chaos in a system of two weakly interacting particles moving in a 
two-dimensional regular confining potential. 

PACS numbers: 05.45.-a, 05.45.Pq, 51.10+y 



1. The chaoticity of motion of the gas of hard scat- 
terers - rigid spheres - was established by N.S. Krylov in 
1943, twenty years before the foundation of the paradigm 
of the modern nonlinear dynamics 1]. He obtained 
the estimate of the Lyapunov exponent a that in two- 
dimensional case has the form 
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where n is the concentration of scatterers, a is the sphere 
(disk) radius, v is some averaged velocity of the particle 
and the condition a^/n <C 1 is fulfilled. This result was 
refined in the following studies of the model of Lorentz 
gas of rigid spheres [2, y| and confirmed numerically j^ . 

More realistic case of soft scatterers, in which the in- 
teraction between the particles smoothly depends on the 
interparticle distance, is much less clear. The result for 
this model derived by Barnett et al. p later was found 
to be in qualitative disagreement with data of numeri- 
cal experiments p. Recently the case of soft scatterers 
has been approached by Kimball Q, but the additional 
superimposed restrictions has led this author just to red- 
erivation of the Krylov estimate Eq. (1). 

In this letter we study the dependence of the Lyapunov 
exponent on parameters of a system of soft scatterers. 

2. We shall use the model of the Lorentz gas - a par- 
ticle of mass m moving with velocity v within a two- 
dimensional array of fixed potential centers located at 
the random points Rj with concentration n. The poten- 
tial energy of the particle at a point r is 



ys(r) = 5]nk-R,l), 
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where the potential of each center will be taken in the 
form V{r) = Vof{ar) with the function f{z) that rapidly 
(e.g. exponentially) decreases with increase of z. In the 



following parameters Vq and a will be called the inter- 
action strength and the inverse interaction range respec- 
tively. 

First, we confine ourselves to the case of low density, 
when the average distance between scattering centers 
L ^ tT^I'^ is much larger than the interaction range: 
hoL ^ 1. Then we can treat the motion of a particle as a 
sequence of collisions with separate centers. The collision 
with the i-th center is characterized by the propagation 
angle 0^ before the collision and the impact parameter of 
the collision pi. 

Next, we assume that the velocity v of the moving 
particle is high enough to justify the use of the small- 
angle approximation, in which the deflection angle Q for 
a collision is given by the expression [^ 
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The linear transformation that relates consequent val- 
ues of small variations of the propagation angle and the 
impact parameter of the following collision can be written 
(see J]) as 
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where Mi is the stability matrix for the i-th collision: 
d{(i>i+i,Pi+i) 
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The form of stability matrix Mi depends essentially on 
position of consequent scattering centers relatively to the 
particle trajectory. This trajectory separates the plane 
in two domains; for convenience we shall refer to them 
as to regions above and below the trajectory. Then we 
have four different cases: 1) both scatterers are above 
the trajectory; 2) the first is above, the second is below; 
3) the first is below, the second is above and 4) both 
scatterers are below the trajectory. For these cases we 



have respectively 
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where ^^ is the derivative of the deflection angle on the 
impact parameter for the i-th collision, 
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and r]i is the distance between the i-th and the (z + 1 
scattering centers. 

The distributions of variables ^i and rji can be charac- 
terized by their first moments, the mean values ^i — (^), 
Vi = iv) s-nd the fluctuations (standard deviations) 
6 = ((C - S.i)Y^^, m = ((77 - 771)^)^/2. The moments of 
^ can be calculated from the distribution of the impact 
parameters of collisions w{p). Since this distribution has 
a characteristic length about the average distance L be- 
tween the scatterers, for L ^ aT^ we can replace this 
distribution by a constant, w(p) -^ k;(0) ~ i-'~^ ■, that 
yields 
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where c\ is a numerical constant of order unity that de- 
pends on the exact form of the interaction potential. The 
distribution of quantities r]i does not depend on the in- 
teraction parameters. Its moments have values 771 — C2L 
and 772 = C3L where C3 and C3 are numerical constants 
of order unity that could be found by means of statisti- 
cal geometry. From the dimensionality considerations it 
follows that the dimensionless Lyapunov exponent "per 
collision" can depend only on dimensionless combinations 
X = |27?i and Y = ^2m- 

Although one can assume the values of £,i and rji to be 
non-correlated, the sequence of types of matrices Mi is 
correlated strongly: all matrices have equal probabilities 
p = 1/4, but after any given one only the matrices of 
two types can follow with probabilities p = 1/2. These 
correlations preclude the use of the standard analytical 
methods of calculation the Lyapunov exponent through 
the small disorder expansion 9] . We note in passing that 
for the ensemble of matrices only of the type M*^^^ or 
only of the type M^'^' with identical rji = fji (772 = 0) 
an analytical expression for the Lyapunov exponent a = 
0.289X^/'^ has been found in the context of the theory of 
one-dimensional Anderson localization [2j. 

Function a{X,Y) has been determined by numerical 
experiments. The sequence of positions of scattering cen- 
ters and values of ^^ (with gaussian distribution) and r/^ 
(with the Poisson distribution) were created by random 
number generators. Then the Lyapunov exponent has 



been calculated for a product of 10^ matrices. Compar- 
ison of different runs shows that this volume of compu- 
tation gives the values of a with an error Act < 5 • 10^*. 
The data found with different parameters (see Fig. 1) 
are accurately fit by the simple linear function 



a = 0.31X 



(9) 
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FIG. 1: Dependence of the Lyapunov exponent a of the 
sequence of random matrices M, given by Eq. (6) on the 
variable X — £,2fji for different values of 771: squares for fji = 
1, filled circles for fji = 2, open circles for fji = 3. The line 
shows the dependence given by Eq. (9). 

Now from Eqs. (8) and (9), taking into account the av- 
erage time interval between collisions Q ~ L/v, we obtain 
for the Lyapunov exponent of a particle moving within a 
two-dimensional array of soft scatterers the expression 
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where c is a numerical constant. This formula presents 
the main result of the paper. The comparison of Eq. 
(10) with the Krylov formula Eq. (1) shows the quali- 
tatilevely opposite dependence on the velocity v and on 
the interaction range a ~ a~^ and quantitatively differ- 
ent dependence on the concentration n. 

3. To check the expression Eq. (10) we consider a 
system of two identical particles moving in a confining 
potential U{x,y). If the motion of a single particle in 
this potential is regular, then chaos in two-particle sys- 
tem can originate only from interaction. If the size A of 
the domain of the configuration space available to parti- 
cles is much larger then the interaction range, Aa 3> 1, 
then we can treat the process of interaction as a series 
of collisions. In this case we expect that the Lyapunov 
exponent will be given by the expression Eq. (10) ap- 
propriately averaged over the distribution of relative ve- 
locities of colliding particles. After this averaging the 
dependence a oc |Vb|\/a must still hold. 

We checked this dependence by a numerical experi- 
ment. The confining potential was taken in the form 
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FIG. 2: Dependence of the Lyapunov exponent cr of a two- 
particle system in a regular confining potential on the inter- 
action strength Vb for the total energy E = 40 and the inverse 
interaction range a = 1. The line shows the fitted function 
o- = 0.0671 Fo I . 
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FIG. 3: Dependence of the Lyapunov exponent cr of a two- 
particle system in a regular confining potential on the inverse 
interaction range a for the total energy E = 40 and the in- 
teraction strength Vo — 1. The line shows the fitted function 
a = 0.042Va. 



U{x,y) — K{x^ -f y^). The regular character of mo- 
tion of a single particle in this potential is secured by 
the existence of two integrals of motion, namely Ix = 
rax^ /2+Kx'^ and ly — my^ /2+Ky'^. The interaction po- 
tential was chosen in the form V{r) = Vo cosh" ar. We 
have used the dimensionless units based on the choice of 
scales m = I, K = 1/4 and the total energy of the system 
E — 40. The dependences of the Lyapunov exponent on 
the interaction strength Vq and on the inverse interaction 
range a are shown in Figs. 2 and 3. 

The accuracy S of agreement between the theory and 
numerical calculations for the dependence a (a) is about 
6%, that is quite satisfactory. However for the depen- 
dence cr(Vo) it is much larger, i5 = 25%, and some expla- 
nation is appropriate. It may be observed that Fig. 2 
demonstrates a persistent asymmetry: numerical values 
of a for the attraction (Vb < 0) happen to be larger (on 
the average by a factor about 1.6) than those for the re- 
pulsion (Vb > 0) with the same |Vb|. In our theory the 
symmetry in sign of Vb comes as a direct consequence of 
the small-angle approximation Eq. (3). From Eq. (10) 
one can infer importance of collisions with small rela- 
tive velocities v, for which the small-angle approximation 
fails. It should be noted that arguments for the asym- 
metry of the observed type were given by Kimball 7] for 
the case of strong scattering. 

In conclusion we note that the model studied in this 
section can serve as a classical counterpart to quantum 
models of two-electron (or electron-hole) systems con- 
fined in quantum dots, that has attracted much attention 
lately 0,Elli3- 
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